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I Dauxois lias recently < 'omim-ul e<l l| an interesting paper by Hilhorst and Schehr [J. 
In the latter paper, two specific probabilistic models, namely a discrete one (from now 
on referred to as the MTG model) and a continuous one 4| (from now on referred to as 
the TMNT model), are analytically discussed. Both models consist in N correlated random 
variables (respectively discrete and continuous), and the point under study is what is the 
limiting distribution when N — > oo. Numerical indications 1 have been found that quite 
strongly suggest that these limiting distributions could be of the (/-Gaussian class (Eq. (1) 
in 1]). Hilhorst and Schehr have shown Q] that they are not, even if they are numerically 
intriguingly close. The relations between the entropic indices q (of the possible g-Gaussian 
limiting distributions) and the parameters of the models have been addressed as well. It 
was analytically confirmed in 2| the correctness of the relations conjectured respectively in 
4|, if these analytically exact non g-Gaussian distributions were to be approached 



3| and 



in 



by g-Gaussians! 

Why may all this have some general interest? The reason lies on the fact that g-Gaussian 

n n □ 

distributions play a special role in nonextensive statistical mechanics [5|, [6|, 17[, a current 
generalization (based on the entropy (2) of [l|]) of the celebrated BG theory. Recent reviews 



of g-statistics and a vast set of experimenta 
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observational and computational applications 
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12]. Very specifically, it has been recently 



and/or verifications can be seen in 

proved Q a g-generalization of the Central Limit Theorem tor </ > 1. More precisely, 
if we have N random variables that are strongly correlated in a special manner (called 
q -independence; see [3| for details), it can be rigorously proved that their sum approaches, 
when appropriately centered and scaled, a g-Gaussian distribution. The standard CLT 
is recovered as the g = 1 particular instance. The proof is under progress for q < 1 
(which are distributions with a compact support, even in the N — > oo limit). Nevertheless, 
preliminary studies suggest that the g-CLT should remain applicable even for g < 1. Since 
the MTG and TMNT models have a compact support (i.e., of the g < 1 type), one expects 
g-Gaussians whenever the correlation is of the g-independent type. Given the results in 



[1] See the title of the present Ref. [3] . 



2j, one is naturally led to argue that the strong correlations in those two models are not 
q-independent but only almost q- independent. An interesting question remains then open. 
What physical ingredient have the MTG and TMNT models failed to incorporate? In other 
words, there is something which is still missing in those models in order for the correlation 
to be exactly q- independent. What is it? Progress is presently being achieved along this 
line (deeply related to asymptotic scale-invariance), but this remains outside the present 
scope, which primarily is the careful scrutiny of the paper 
the weaknesses of that paper. 



l|. Let us address now some of 



(i) In contrast with what one reads in 1[, the TMNT model is not at all described 
nor presented in its Ref. [3]. In fact, it is not even mentioned there, and has never been 
published! 2 . Dauxois visibly confuses with [14j, whose content has absolutely nothing to 
do with the model discussed by Hilhorst and Schehr. Had the Author of [1] paid more 
attention to his own Ref. [3], or at least to its title, the confusion would have been avoided. 



(ii) The statement "could be the basis for a generalized central limit theorem" (below Eq. 
JJI) reveals that the Author is possibly unaware of the proof existing since already 

n nn 

some time in [13j (further extensions are possible along the lines of [la, Il6l]). 



(1) of 



(iii) The statement "The basis for this suggestion ..." (above Eq. (2) of [l|) reveals 
that the Author is unaware that the crucial point is not the isolated fact that o-Gaussians 

n 

optimize S q (something which, in contrast with what is referred in [l|, is not even vaguely 



mentioned in Ref. [21 of 

nn 

to be [17|, [18[] exact stab 
decades called Porous Medium Equation). The mere fact that g-Gaussians would optimize 



but rather the remarkable fact that g-Gaussians also happen 
to be [ill [lsj] exact stable solutions of the nonlinear Fokker-Planck equation (since many 



[2] The details of the TMNT model, as well as all our numerical and graphical results, were transmitted 
privately by me to Hilhorst, who had heard, in Natal-Brazil in March 2007, an oral presentation of mine 
including that subject 
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some specific entropic functional would certainly not be of any particular significance. The 
entire rationale of what was at the time a conjecture (and is now a proved theorem) can be 
seen in 19]. 



of 



(iy) The qualification "in the absence of firm grounds..." reveals that the Author 
lj is possibly unaware of the numerous rigorous results which precisely provide a 



firm mathematical basis for the entropy S„ and the associated nonextensive statistical 



mechanics. Among many others, let us mention the g-generalization, respectively in 



20) 



and in 2l| . of Shannon's and of Khinchine's uniqueness theorems, the already mentioned 



q-central limit theorems in 



13| . the Lesche-stabi 



231 ] of the entropy S q , the analytical results in 24J for Langevin-like stochastic equations 



ity 



221 ] and the Topsoe-factorizability 
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2fl 



in the presence of multiplicative noise, the analytical results of 
unimodal dissipative one-dimensional maps, the analytical connection [29 
Beck-Cohen superstatistics, etc. 
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28 
30, 



concerning 



3l| to the 



(v) The Author of l] wonders whether the whole idea of g-statistics could be "just a 
nice idea and a powerful fitting function". Fitting functions typically have one or more 
fitting parameters. Let us focus on this point. Fitting parameters can be of quite different 
natures. The first, and in some sense more fundamental, class is constituted by those 
numbers for which no theory exists in contemporary physics that could produce them 
from more basic principles. In this class we have all the so-called universal constants of 
contemporary physics (c, h, G, and which, as well known, might not even be constants! 
They are obtained, nowadays with great precision, from fittings concerning basic physical 
laws such as the Planck law for the black-body radiation. A second class is constituted by 
those parameters that can ideally be deduced from first principles through some specific 
way involving some specific information. But, in many cases, the procedure turns out to be 
untractable, either because of mathematical difficulties, or because some of the necessary 
information is missing. Such is the case of the orbit of say Mars. Newtonian mechanics 
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can in principle predict this orbit given, at some time, the positions and velocities of all 
masses of the planetary system. Since this information is unavailable (and even if it was 
available, a colossal computer would be needed to perform the calculations!), astronomers 
use, as a first approximation, the elliptic form of a Keplerian orbit and just fit this form 
(which can be easily obtained within Newtonian mechanics) in order to obtain the concrete 
values of the two axis of the ellipse. The entropic index q belongs to this class. It can be 
in principle deduced from the microscopic dynamics of the system. And this is indeed done 
successfully in some few simple cases. But in most cases, this precise dynamics is unknown, 
or the calculations are untractable. Therefore q can be obtained by fitting results with the 
analytic forms (g-exponentials, g-Gaussians) that emerge within nonextensive statistical 
mechanics. A third class of fitting parameters would be those which emerge from no specific 
theory, but are present in some convenient heuristic expressions. Within the ^-statistical 
theory, we are not especially interested in those. 



(vi) Dauxois considers "particularly pressing" an undoubtedly relevant question, namely 
"does the q-Gaussian law describe the details of some physical problems?" . Many systems 
exist that have been shown to exhibit, at the available precision, g-Gaussians. A particularly 
instructive example concerns cold atoms in dissipative optical lattices. It was predicted in 



2003 by Lutz 32] that these atoms should have a g-Gaussian distribution of velocities with 



q 



(i) 



Er and Uq being respectively the microscopic recoil energy and potential depth. In 2006, an 
impressive confirmation was provided 



331 ] , through both quantum Monte Carlo simulations, 
and real experiments with Cs atoms (the quality of the agreements are given by the 
correlation coefficients R 2 = 0.995 and R 2 = 0.9985 respectively). Other examples (some of 
them briefly mentioned in l|) that, within reasonably good approximation, have exhibited 



g-Gaussians are the motion of cells of Hydra viridissima 
fluctuations of the magnetic field in the solar wind 



temperature of the universe 
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38, 



39 



42], financial return distributions 



34 1, defect turbulence |35j, 



40 



44 



, fluctuations of the 



45 



46, 



47|, 



48 



49|, 



many-body Hamiltonian with long-range-interacting classical rotators approached through 
the HMF [50, |5l| , logistic map at edge of chaos [52j , granular matter [53j , earthquakes |54| , 
etc. 



(vii) "More importantly", as a second, "particularly pressing" question, Dauxois asks "is 
anyone able to provide analytical predictions of the value of the q-index in terms of the 
microscopic parameters of the physical system?". This is of course an important point, that 
any quantitatively mature statistical mechanical theory (and physical theory in general) is 
expected to satisfy. The Author might be unaware of the fact that many examples exist 
for which analytical calculations of q have been possible from microscopic or mesoscopic 
dynamics. Let us mention some of them. Of course, Eq. (1) already constitutes one such 



55) 



example. Another analytical result is the fact that the index q has been connected 
with the standard critical exponent 6 (the exponent characterizing the order parameter as a 
function of its thermodynamically conjugate variable, at a second-order critical phenomenon) 
as follows: 



In the area of quantum entanglement, the block entropy of fully entangled quantum ma; 



chains involving only short-range interactions, has been recently proved to satisfy 



gnetic 



56] 



, = , (3) 

where c is the central charge emerging within conformal quantum field theory. This means 
that, for the Ising ferromagnet (i.e., c = 1/2) we have q = y/Wf — 6 = 0.08..., and for the 
isotropic XY ferromagnet (i.e., c = 1), we have q = vTO — 3 = 0.16.... In the limit c — ■> oo 
the BG result q = 1 is recovered. As another analytical connection, we may mention the 



Albert-Barabasi model in 



see 



for (asymptotically) scale-free networks, for which we have 



58j for the connection of the Albert-Barabasi exponent and q) 



2m(2 — r) + 1 — p — r 

Q = —To o A i 1 ' 4 ) 

m(o — Ir) + 1 — p — r 
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where fm, p, r) are microscopic parameters of the model. Let us also mention that, on the 



grounds of 



18| . the following connection is since long predicted between q and the anomalous 



diffusion space-time scaling exponent 7: 



7 



q 



(5) 



This prediction has been verified by now in various systems, such as the Hydra viridissima 



cells |34|, defect turbulence |35j, a-XY inertial ferromagnetic model [59J (see [60| for details 



on the model), and granular matter flowing out from a silo 



53] 



As our last present illustration of analytical predictions involving q, we may mention one 



which is available since 1999 



631 ] and which has been successfull y appli ed to plasma 



astrophysics (e.g., solar plasma). For weakly nonideal plasma, we have 



61 



62 



q-l\ = 24a 4 7 2 < 1 



(6) 



where a (typically 0.4 < a < 0.9) is the ion-ion correlation-function parameter (see details 
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65(|). and 7 = (Ze) 2 n 1//3 /ksT is the plasma parameter (see Eq. (7) of 6l|), n being 



the average density. 



(viii) The Author of lj makes, a few lines later, a short-cut which might mislead the 
quick reader. He states that "Hilhorst and Schehr show that q-Gaussians do not pass 
a careful inspection. ". Instead, what these Authors have (interestingly) shown is that 
the MTG and TMNT models do not pass in the very specific sense that their limiting 
distributions are not exactly but only approximatively g-Gaussians. This point was already 
discussed in the beginning of the present manuscript. 



fix) The Author of [l| states that "if there is one lesson that has to be learned here, 
it is that one should be extremely careful when interpreting non- Gaussian data in terms 
of q-Gaussians". This is certainly true. But, of course, this also is true for Gaussians, 
Lorentzians, circles, ellipses, straight lines (in fact, g-Gaussians are mere straight lines 
in hi g y = (y l ~ q — 1)/(1 — g) versus x 2 representation), parabolas, etc, all the analytical 



tools with which physicists have attempted to interpret and quantitatively understand the 
plethora of experimental, observational and computational results they are dealing with 
since centuries and millenia. But of course, - - ce qui est clair sans le dire est encore 
plus clair en le disant! — , any experimental, observational or computational setup can 
only provide a finite amount of numbers with finite precision. This will always allow for 
an infinite number of analytic curves or interpretations. The content of [2j provides an 
interesting, and possibly fruitful, illustration of this obvious fact. 



Acknowledgements 

Financial support by CNPq and Faperj (Brazilian Agencies) is gratefully acknowledged. 



[1] T. Dauxois, Non-Gaussian distributions under scrutiny, 0712.3153 [cond-mat.stat-mec]; J. 
Stat. Mech. (2007) N08001. 

[2] H.J. Hilhorst and G. Schehr, A note on q-Gaussians and non-Gaussians in statistical mechan- 
ics, J. Stat. Mech. (2007) P06003. 

[3] L.G. Moyano, C. Tsallis and M. Gell-Mann, Numerical indications of a q-generalised central 
limit theorem, Europhys. Lett. 73, 813 (2006). 

[4] W. Thistleton, J.A. Marsh, K. Nelson and C. Tsallis, unpublished (2006). 

[5] C. Tsallis, Possible generalization of Boltzmann-Gibbs statistics, J. Stat. Phys. 52, 479 (1988). 

[6] E.M.F. Curado and C. Tsallis, Generalized statistical mechanics: connection with thermody- 
namics, J. Phys. A 24, L69 (1991); Corrigenda: 24, 3187 (1991) and 25, 1019 (1992). 

[7] C. Tsallis, R.S. Mendes and A.R. Plastino, The role of constraints within generalized nonex- 
tensive statistics, Physica A 261, 534 (1998). 

[8] M. Gell-Mann and C. Tsallis, eds., Nonextensive Entropy - Interdisciplinary Applications, 



8 



(Oxford University Press, New York, 2004). 
[9] J. P. Boon and C. Tsallis, eds., Nonextensive Statistical Mechanics: New Trends, New Perspec- 
tives, Europhysics News 36 (6) (European Physical Society, 2005); Corrigendum: Europhysics 
News 37 (1), 25 (2006). 

[10] S. Abe, H.J. Herrmann, P. Quarati, A. Rapisarda and C. Tsallis, eds. Complexity, Metastability 
and Nonextensivity, American Institute of Physics Conference Proceedings 965 (New York, 
2007). 

[11] C. Tsallis, Entropy, in Encyclopedia of Complexity and Systems Science (10 volumes) 

(Springer, Berlin, 2007), in press. 
[12] A regularly updated bibliography is available at http://tsallis.cat.cbpf.br/biblio.htm 



[13] S. Umarov, C. Tsallis and S. Steinberg, A generalization of the central limit theorem consis- 
tent with nonextensive statistical mechanics, cond-mat/0603593 (2006); C. Tsallis and S.M.D. 
Queiros, Nonextensive statistical mechanics and central limit theorems I - Convolution of inde- 
pendent random variables and q-product, in Complexity, Metastability and Nonextensivity, eds. 
S. Abe, H.J. Herrmann, P. Quarati, A. Rapisarda and C. Tsallis, American Institute of Physics 
Conference Proceedings 965, 8 (New York, 2007); S.M.D. Queiros and C. Tsallis, Nonexten- 
sive statistical mechanics and central limit theorems II - Convolution of q-independent random 
variables, in Complexity, Metastability and Nonextensivity, eds. S. Abe, H.J. Herrmann, P. 
Quarati, A. Rapisarda and C. Tsallis, American Institute of Physics Conference Proceedings 
965, 21 (New York, 2007). 

[14] W. Thistleton, J. A. Marsh, K. Nelson and C. Tsallis, Generalized Box-Muller method for 
generating q-Gaussian random deviates, IEEE Transactions on Information Theory 53, 4805 
(2007). 

[15] S. Umarov, C. Tsallis, M. Gell-Mann and S. Steinberg, q- generalization of symmetric a-stable 



distributions. Part I, |cond-mat/ 0606038 (2006); S. Umarov, C. Tsallis, M. Gell-Mann and 
S. Steinberg, q- generalization of symmetric a-stable distributions. Part II, cond- mat/0606040 
(2206). 



9 



[16] S. Umarov and C. Tsallis, On multivariate generalizations of the q-central limit theorem consis- 
tent with nonextensive statistical mechanics, in Complexity, Metastability and Nonextensivity, 
eds. S. Abe, H.J. Herrmann, P. Quarati, A. Rapisarda and C. Tsallis, American Institute of 
Physics Conference Proceedings 965, 34 (New York, 2007); C. Vignat and A. Plastino, Central 
limit theorem and deformed exponentials, J. Phys. A 40, F969 (2007). 

[17] A.R. Plastino and A. Plastino, Non-extensive statistical mechanics and generalized Fokker- 
Planck equation, Physica A 222, 347 (1995). 

[18] C. Tsallis and D.J. Bukman, Anomalous diffusion in the presence of external forces: exact 
time- dependent solutions and their thermostatistical basis, Phys. Rev. E 54, R2197 (1996). 

[19] C. Tsallis, Nonextensive statistical mechanics, anomalous diffusion and central limit theorems, 
Milan Journal of Mathematics 73, 145 (2005). 

[20] R.J.V. dos Santos, Generalization of Shannon's theorem for Tsallis entropy, J. Math. Phys. 
38, 4104 (1997). 

[21] S. Abe, Axioms and uniqueness theorem for Tsallis entropy, Phys. Lett. A 271, 74 (2000). 

[22] S. Abe, Stability of Tsallis entropy and instabilities of Renyi and normalized Tsallis entropies: 
A basis for q-exponential distributions, Phys. Rev. E 66, 046134 (2002). 

[23] F. Topsoe, Factorization and escorting in the game-theoretical approach to non-extensive en- 
tropy measures, Physica A 365, 91 (2006). 

[24] C. Anteneodo and C. Tsallis, Multiplicative noise: A mechanism leading to nonextensive 
statistical mechanics, J. Math. Phys. 44, 5194 (2003). 

[25] F. Baldovin and A. Robledo, Sensitivity to initial conditions at bifurcations in one- dimensional 
nonlinear maps: Rigorous nonextensive solutions, Europhys. Lett. 60, 518 (2002). 

[26] F. Baldovin and A. Robledo, Universal renormalization- group dynamics at the onset of chaos 
in logistic maps and nonextensive statistical mechanics, Phys. Rev. E 66, R045104 (2002). 

[27] F. Baldovin and A. Robledo, Nonextensive Pesin identity. Exact renormalization group an- 
alytical results for the dynamics at the edge of chaos of the logistic map., Phys. Rev. E 69, 
045202(R) (2004). 



10 



[28] E. Mayoral and A. Robledo, Tsallis' q index and Mori's q phase transitions at edge of chaos, 
Phys. Rev. E 72, 026209 (2005). 

[29] G. Wilk and Z. Wlodarczyk, Interpretation of the nonextensivity parameter q in some appli- 
cations of Tsallis statistics and Levy distributions, Phys. Rev. Lett. 84, 2770 (2000). 

[30] C. Beck, Dynamical foundations of nonextensive statistical mechanics, Phys. Rev. Lett. 87, 
180601 (2001). 

[31] C. Beck and E.G.D. Cohen, Superstatistics, Physica A 322, 267 (2003). 

[32] E. Lutz, Anomalous diffusion and Tsallis statistics in an optical lattice, Phys. Rev. A 67, 
051402(R) (2003). 

[33] P. Douglas, S. Bergamini and F. Renzoni, Tunable Tsallis distributions in dissipative optical 
lattices, Phys. Rev. Lett. 96, 110601 (2006). 

[34] A. Upadhyaya, J.-P. Rieu, J. A. Glazier and Y. Sawada, Anomalous diffusion and non-Gaussian 
velocity distribution of Hydra cells in cellular aggregates, Physica A 293, 549 (2001). 

[35] K.E. Daniels, C. Beck and E. Bodenschatz, Defect turbulence and generalized statistical me- 
chanics, Physica D 193, 208 (2004). 

[36] L.F. Burlaga and A.F. Vinas, Multiscale structure of the magnetic field and speed at 1 AU 
during the declining phase of solar cycle 23 described by a generalized Tsallis PDF, J. Geophys. 
Res. - Space Phys. 109, A12107 (2004). 

[37] A.F. Vinas, R.L. Mace and R.F. Benson, Dispersion characteristics for plasma resonances 
of Maxwellian and Kappa distribution plasmas and their comparisons to the IMAGE/RPI 
observations, J. Geophys. Research 110, A06202 (2005). 

[38] L.F. Burlaga and A.F.-Vinas, Triangle for the entropic index q of non-extensive statistical 
mechanics observed by Voyager 1 in the distant heliosphere, Physica A 356, 375 (2005). 

[39] L.F. Burlaga, N.F. Ness and M.H. Acuna, Multiscale structure of magnetic fields in the he- 
liosheath, J. Geophysical Res.-Space Physics 111, A09112 (2006). 

[40] L.F. Burlaga, A. F- Vinas and C. Wang, Tsallis distributions of magnetic field strength varia- 
tions in the heliosphere: 5 to 90 AU, J. Geophysical Res.-Space Physics 112, A07206 (2007). 



11 



[41] L.F. Burlaga, N.F. Ness, M.H. Acuna, Magnetic fields in the heliosheath and distant helio- 

sphere: Voyager 1 and 2 observations during 2005 and 2006, Astrophys. J. 668, 1246 (2007). 
[42] A. Bernui, C. Tsallis and T. Villela, Deviation from Gaussianity in the cosmic microwave 

background temperature fluctuations, Europhys. Lett. 78, 19001 (2007). 
[43] L. Borland, Closed form option pricing formulas based on a non-Gaussian stock price model 

with statistical feedback, Phys. Rev. Lett. 89, 098701 (2002). 
[44] L. Borland, A theory of non-gaussian option pricing, Quantitative Finance 2, 415 (2002). 
[45] S.M.D. Queiros, On the emergence of a generalised Gamma distribution. Application to traded 

volume in financial markets, Europhys. Lett. 71, 339 (2005). 
[46] S.M.D. Queiros, On non- Gaussianity and dependence in financial in time series: A nonexten- 

sive approach, Quant. Finance 5, 475-487 (2005). 
[47] J. de Souza, L.G. Moyano and S.M.D. Queiros, On statistical properties of traded volume in 

financial markets, Eur. Phys. J. B 50, 165-168 (2006). 
[48] S.M.D. Queiros, L.G. Moyano, J. de Souza and C. Tsallis, A nonextensive approach to the 

dynamics of financial observables, Eur. Phys. J. B 55, 161 (2007). 
[49] S.M.D. Queiros, Are all highly liquid securities within the same class?, Eur. Phys. J. B 60, 

265 (2007). 

[50] A. Pluchino, A. Rapisarda and C. Tsallis, Nonergodicity and central limit behavior in long- 
range Hamiltonians, Europhys. Lett. 80, 26002 (2007). 

[51] C. Tsallis, A. Rapisarda, A. Pluchino and E.P. Borges, On the non-Boltzmannian nature of 
quasi- stationary states in long-range interacting systems, Physica A 381, 143 (2007). 

[52] U. Tirnakli, C. Beck and C. Tsallis, Central limit behavior of deterministic dynamical systems, 
Phys. Rev. E 75, 040106(R) (2007). 

[53] R. Arevalo, A. Garcimartin and D. Maza, Anomalous diffusion in silo drainage, Eur. Phys. J. 
E 23, 191 (2007). 

[54] F. Caruso, A. Pluchino, V. Latora, S. Vinciguerra and A. Rapisarda, Analysis of self- organized 
criticality in the Olami-Feder-Christensen model and in real earthquakes, Phys. Rev. E 75, 



12 



055101(R)(2007). 

[55] A. Robledo, Unorthodox properties of critical clusters, Molecular Physics 103, 3025 (2005). 

[56] F. Caruso and C. Tsallis, Extensive nonadditive entropy in quantum spin chains, in Complex- 
ity, Metastability and Nonextensivity, eds. S. Abe, H.J. Herrmann, P. Quarati, A. Rapisarda 
and C. Tsallis, American Institute of Physics Conference Proceedings 965, 51 (New York, 
2007). 

[57] R. Albert and A.-L. Barabasi, Topology of evolving networks: Local events and universality, 
Phys. Rev. Lett. 85, 5234 (2000). 

[58] S. Thurner, F. Kyriakopoulos and C. Tsallis, Unified model for network dynamics exhibiting 
nonextensive statistics, Phys. Rev. E 76, 036111 (2007). 

[59] A. Rapisarda and A. Pluchino, Nonextensive thermodynamics and glassy behaviour in Hamil- 
tonian systems, Europhysics News 36, 202 (2005) [Erratum: 37, 25 (2006)]. 

[60] C. Anteneodo and C. Tsallis, Breakdown of the exponential sensitivity to the initial conditions: 
Role of the range of the interaction, Phys. Rev. Lett. 80, 5313 (1998). 

[61] M. Coraddu, G. Kaniadakis, A. Lavagno, M. Lissia, G. Mezzorani and P. Quarati, Thermal 
distributions in stellar plasmas, nuclear reactions and solar neutrinos, in Nonextensive Statis- 
tical Mechanics and Thermodynamics, eds. S.R.A. Salinas and C. Tsallis, Braz. J. Phys. 29, 
153 (1999). 

[62] M. Coraddu, M. Lissia, G. Mezzorani and P. Quarati, Deuterium burning in Jupiter interior, 

Physica A 305, 282 (2002). 
[63] M. Lissia and P. Quarati, Nuclear astrophysical plasmas: Ion distributions and fusion rates, 

Europhysics News 36, 211 (2005). 
[64] X.Z. Yan and S. Ichimaru, Phys. Rev. A 34, 2167 (1986). 

[65] S. Ichimaru, Statistical Plasma Physics, Vol. 1, Paragr. 1.3 (Addison Wesley, Redwood City, 
1992). 



13 



